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1 Introduction 

^ ■ . . 

The inverse scattering transform method (ISM) is an effective tool for study- 
^ . ing nonhnear integrable equations [[H . It allows one to construct large classes 



^ ; of exact solutions and to investigate in detail the Cauchy problem includ- 

^ ing the writing down explicit time asymptotics. As for the mixed problem 

> ! for the same equations, here the success of the method is not so impressive. 

O ; For instance, the ISM applied to the initial boundary value problem in its 

^ ; general formulation is not sufficiently effective (see., for example, 0), except 

^ '■ some special kinds of boundary conditions (see. 0, 1, 0, The problem of 

searching boundary conditions consistent with the integrability property of 
the equation given as well as the problem of finding proper modifications of 
O ■ analytical integration procedure for the corresponding boundary value prob- 

lems is undoubtedly important. 
^ ; Notice that in all known cases of initial boundary value problems studied 

up by means of ISM the integration algorithm is based on reduction to the 
Cauchy problem. One prolongs the initial value by means of the Backlund 
transform's ordinary differential equations. This way is available if the equa- 
tion is invariant under a reflection type transformation x ^ —x only. 

In our articles [jTj], a symmetry test was worked out to verify whether 
the boundary value problem is consistent with integrability or not. Utilizing 
this test in [g] and [Q] some examples of boundary conditions were found 
for the Korteweg-de Vries type equations. Evidently these equations don't 
admit any x reflection. So to integrate the corresponding initial boundary 
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value problems the continuation algorithm mentioned above is not suitable 
and in the present paper an alternative approach is proposed. The essence 
of our approach consists of complementing the conjugation contour of usual 
associated Riemann problem by adding new branches. As an illustrative 
example we consider the following initial boundary value problem for the 
Modified KdV equation 

ut = Uj;xx + Geu^u^, £ = ±1, (1) 

u\or^O = 0, Ua:\x^O = 0, (2) 
u\t=Q = Uo{x), Uo{x)\j:^+oo ^ (3) 

in the quadrant x > 0, t > 0. The consistency condition uq{0) = woa;(0) = 
of the initial value with the boundary conditions is supposed to be valid at 
the corner point a: = 0, t = 0, and besides the initial value uo{x) is required 
to vanish fast enough. The problem under consideration is well posed. Notice 
that the correctness of mixed problems like (|I])-(0) for linear and quasi-linear 
third order equations is a well studied question (see, for instance, [|TD|)- It is 



worthwhile to say that the similar problem on the left semi-line a: < is not 
correctly posed - in this case one should put only one boundary condition at 
the point x = 0. This fact implies immediately that the problem ([I|)-(0) can 
not be reduced to the Cauchy problem on a real axis when the initial value 
uo{x) is chosen arbitrary. 

Additionally we put rather severe requirement upon the initial data mean- 
ing the absence of discrete eigenvalues of the auxiliary linear system of equa- 
tions (see conjectures i)-iv) in section 5 below) thus keeping out of our at- 
tention an important question on soliton-like solutions to the boundary value 
problem (|T])-(0)- In this connection mention the recent article []TT|, were an 
explicit description of algebro-geometric solutions to the KdV equation is 
given satisfying integrable boundary conditions. 

Discuss briefly the content of the paper. In the second section we consider 
a linear analogue of the problem (|I])-(0) and construct its explicit solution by 
the Wiener-Hopf method. In third one the Lax representation is formulated 



for the initial boundary value problem on a half-line. In fourth - the direct 
scattering problem is considered for the associated linear system of equations. 
On a half-line solution of this problem is not unique in contrast to the case 
of the whole line []T^. Some freedom arises even in choosing the conjugation 
contour - here it does not coincide with the continuous spectra of the x- 
equation. In §5 it is proved that the scattering matrix defines uniquely the 
coefficients of equation. The inverse scattering problem is reduced to the 
Riemann problem posed on six rays dividing the complex plane to six equal 
sectors. In §6 the time evolution of scattering data is defined and it is shown 
that potential recovered in inverse problem solves the initial boundary value 
problem ([I|)-(^. In the last seventh section the solvability of the Riemann 
problem on system of rays is discussed. 



2 Solution of mixed problem in linear case 

This section is an auxiliary one, its aim is to explain with a simple illustrative 
example the main parts of the algorithm proposed and in such a way make 
easier the reading of the next sections. 

Consider a linear version of the problem ([I|)-(0): 

Ut = Uxxx, (4) 
U\x={) = 0, Ux\x^O = 0, (5) 
u\t=o = qix), 9(a:)|^^+oo ^ 0. (6) 

Apply to the equation (g) the Fourier transform instead of Laplace one as it 
used to be. Set v{^^t) = Jq° u{x,t) exp{i^x)dx, then apparently the function 
v{^,t) solves the equation 

vt{tt)=u,,{0,t) + iev{^,t), (7) 

which contains an unknown function Uxxi^^t)- In order to exclude the un- 
certainty arisen let us do a change of variables in the equation putting 

2,i7r . 

h{^,t) = v{^,t) — v{u^,t), where u = e~ is a cubic root of the unity. The 
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function h{^, t) satisfies much more simple equation ht = i^^h which is solved 
trivially. Now it is necessary to invert the change of variables. To this end as 
well as in order to solve the scattering problem below we need in the follow- 
ing geometric figures. Define on the complex plane six congruent sectors /i, 
/2, ... /e, by dividing the complex plane by six rays ended at the origin: 
^j = {i^, CLrg{^) = |(j — 1)} enumerated in such a way that Ij is a sector 
located between rays ij ij+i (see picture below) 




Picture. Conjugation contour. 

Note that by construction the function h{^,t) for t = is analytical in sector 
Ii and vanishes at infinity. These properties are valid also for i > 0, but 
violated for t < 0. One can check directly that the function searched is 
represented for t > as a sum of three Cauchy type integrals 

2™K) = f ^ + f + f (8) 

Jii S — C, "'*4 s — ^ s — t, 

Here the integration is taken along rays £j in direction from the origin to the 
infinity. The condition /i(0) = guarantees the continuity of the function (|8D 
in the closed upper half-plane. 
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The formula is simplified essentially when ^ G /i U /2. Really for such 
^ integration along the ray £4 can be replaced by integration along £3, and 
afterwards by means of the change of variables i' = one can reduce sum of 
three integrals in to a integral along £1: 

2mv{^)= ^ 

1 z/ — 4 

For ^ G /s the function v{^) is found from: = v{ix}'^^) — h{uP'^). Applying 
the inverse Fourier transform to v{^^ t) one finds the solution searched t) 
to the problem (g)-(J5|). 



3 Lax representation for the problem on a half-line 

Consider for convenience a system of equations in partial derivatives contain- 
ing equation ([T]) as a particular case 

Ut = Ur^xx ~ QuxUu'^, ul = ul^.^^ — Gu^uu^. (9) 

Put in accordance with (g) and (g): 

U\x^0 = Ux\x^Q = U'^lx^Q = U^lx^Q = (10) 

and 

u{x,0) = uq{x) ^ 0, u^{x^O) = Uq{x) ^ for x ^ 00. (11) 

All our constructions below are consistent with the involution = —su re- 
ducing the problem (|i)-(|n|) to the problem (|I])-(§). The system of equations 



is a compatibility condition of the following over determined system of 
equations 

yx = H^<J3 + U)y, (12) 
yt = m'<J3 - - ^V, - Vo)y, (13) 



where 



-2iuu^ —2u 
— 2w+ 2mw" 



x 
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iu \ 






j 




V — m+ 



Along the line a: = by means of the boundary conditions (|rD|) equation ([TBj) 
takes the form 

yt = (4^^V3 - Uxx)y. (14) 

Changing parameter u = — 4^*^ in the last equation one gets again Zakharov- 
Shabat equation of the form ([T2|). 

Thus three linear systems of equations ([T^) , ([T^) and provide the Lax 
representation for the initial boundary value problem (0), ([TD|) and ([TT]) . 



4 Direct scattering on a half-line 

Using the fact that the potential U{x^O) vanishes rapidly enough one can 
construct a matrix valued solution y{x,^) of the auxiliary system ([12]), by 
prescribing for x +oo the following asymptotically trigonometric behavior 

valid for any real ^. As it is known such a matrix valued solution always 
exists and its columns y^^''^\x, ^) define two vector valued functions il){x, ^) = 
y^{x,^)e'^^^ and = ^)e^*^^, admitting analytic continuations with 

respect to the parameter ^ from the real axis onto the lower and, respectively, 
upper half-planes. Moreover, the matrix valued function 

tends to the unit matrix when R 3 ^ ^ oo. The matrix T{^) = ('?/^(^), 0(0)' 
compound by columns = i/j{0^^) and = 0(0,0 is called scattering 
matrix of the system ([T2| ) on a half line. 

Since to solve completely the direct scattering problem means to find a 
pair of linearly independent piece- wise analytic vector valued solutions to the 
system of equations (fT^l) which should be defined at all points of the ^-plane 
except may be a set of measure equal to zero (see, for instance, []T^), hence we 
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need in one more solution. In our case the variable x takes only non-negative 
values, so we have a freedom in choosing the second analytic solution. Utilize 
this freedom in the following way. Compound matrix valued functions ^'j(^), 
analytic in sectors 7^, with the help of columns of the scattering matrix T{^) 
by putting dawn 

Complete the discussion on the direct scattering problem defining the solution 
y{x, of the linear problem (|r2|), by giving up its value at the point x = 0, 

z/(o,0 = *.(0, ^^ij- 

One can prove that the function ^(x,^) = given as 

is a piece-wise analytic function, defined on the whole complex ^-plane. In 
passing up from a sector Ij to the sector Ij+i the ?/;-function have a jupm: 

= ^,(2:, 0^.(^,0- (16) 

The system of correlations (fTOj ) can be considered as the Riemann problem 
on rays {ij}. Solutions to the problem ( [T^ ) satisfy some additional constraint 
posed at the infinity oo) = 1. Besides the following asymptotic repre- 

sentation can be shown to be valid 

N 

k=i 

if the coefficients uo{x), Uq{x) of system of equations (|T2]) are smooth, rapidly 
decreasing functions. It can easily be checked that the conjugation matrices 
depend on the variable x in an explicit way 

Rj{x,0=e-'^^^'RMOe^^^^\ (18) 
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5 Inverse scattering problem 



In the inverse scattering the problem is considered to recover the potential 
U{x) of the system (|T2|) for a scattering matrix given. Remind (see above) 
that the scattering matrix is described by a pair of vector functions '?/^(^), 0(0? 
analytical in lower and, respectively, upper half-planes. Coordinates of these 
vectors are connected by equation V^i(^)02(O ~ ^2(0'5^i(0 = 1 for ^ G R. 
Because of the analytical properties of entries this equation allows one to 
express functions V^i(O) 02(0 terms of V^2(0 0i(O by well known 
formulae (see if zeros of functions 02 (0? located in domains Im^ < 

and Im^ > 0, respectively, are hsted a priori, and 1 + V^2(O0i(O for 
Im^ = 0. 

Note that the scattering matrix of the half-line problem defines uniquely 
the potential of system of equations ([T^) in contrast to the case of whole line. 
Let us formulate a precise statement. 

Proposition 1. Let T{^) and T{^) are scattering matrices, correspond- 
ing the potentials U{x) U{x). Then the equation T{^) = T{^) implies the 
equation U{x) = U{x). 

Scheme of proof. Suppose that T{^) = T{^). It means that the Riemann 
problem (|TB|) has two different solutions ^j+i(x,0 = '^j{x,^)Rj{x,^) and 
^j+i{x,^) = '^j{x,^)Rj{x,^). Let f{x,^) be a ratio of these two solutions 

/(x,o = ^',(x,o*7'(:c,0 = ^',+i(x,0*7+i(^,0- 

Evidently function f{x,^) is bounded for large values of ^, satisfies a linear 
equation with respect to x: fx = i^[f, cr^] + ufa2 — ua2f, besides /(O, ^) = 1. 
Consequently, it is analytical on whole complex plane and according to the 
Liouville theorem it doesn't depend on ^. It follows from the linear equation 
on / that [/, CTs] = 0. Thence, fx = and therefore f{x,^) = 1, that means 
that U{x) = U{x). 

Exclude from our consideration some degenerate cases in the inverse scat- 
tering problem by imposing the following additional requirements: 

l + V^2(O0i(O7^O for ^gR, 



n) l0g(l + ^2(001 (0)11=1^ = 0, 

III) (0(O,V^(a;'O)7^O for ^G/iU/2, 
vi) (0(0, ^K)) 7^0 for ^G/sUJa, 



where the skew-scalar product (X, Y) between two-dimensional vectors X 
and Y is understood as the determinant of matrix: (X, Y) = det(X, Y) and 
the bar over a letter denotes closure of a multitude. 

Proposition 2. Assume we are given a pair of functions 0i(O? "02(0 such 
that 



where fi^2{x) vanish fast enough together with all derivatives. Recover the 
matrix T{^) = (i/ji^), 0(0)' choosing the first coordinate of the vector i/j and 
the second coordinate of the vector as Cauchy type integrals 



Let the all conditions i)-vi) are valid. Then matrix ^(0 is a scattering matrix 
of the problem on a half-line if and only if for all a: > a non-degenerate 
solution (det ^ for all ^ G /j) exists to the Riemann problem ([T^) with 
the conjugation matrices defined as 









i?2 = -R4 = 1 
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where 



Scheme of proof. The necessity of the conditions ([T^), ([T^) was proven 



in the previous section. Prove now the sufficiency. Define the conjugation ma- 
trices Rj for the scattering matrix T{^) = ("0(0, 0(0) given as it is pointed 
out in the proposition. Suppose that for aU x > the Riemann problem ([TBI), 
(|T8D is uniquely solvable. Find the potential U (x) of the system (|T2|). To this 



end differentiate equation (|T6|) with respect to x and then transform it to the 
form 

Applying now the Liouville theorem and asymptotic representation (\n\) it is 
easy to get that g = —i^ds + t/, i.e. 

= -^0^3, + (19) 

Remark to the proposition 2. The functions Aij and conjugation 
matrices Rj admit some analytical properties. For instance, An is analytic 
in the domain /i U /e, A12 - in the sector /i, ^421 - in the sector I2 and, 
lastly, ^422 is analytic in /1U/2. Conjugation matrices Ri and R5 are analytic, 
respectively, in I2 and Ik,. The matrices R^ and Rq are defined, generally 
speaking, only on corresponding contours £4 and £1. 



6 Time dynamics and boundary conditions 

The scattering matrix evaluates in time in a very complicated way - it is the 
main difficulty in investigating the initial boundary value problem. One can 
get some impressions about the problem by looking at the equation (|7|) on the 
function v{^, t) in the linear case. On the other hand side the conjugation ma- 
trices which are some rational functions of the elements of the scattering 
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matrix, depend on t explicitly. It is more convenient to postulate time dynam- 
ics in the inverse scattering problem than to find it up in the direct problem. 
Let the conjugation matrices Rj{x,^,t) of the Riemann problem (|TB|), (fTSI) 
evaluate in time by means of the following equations Rj^t = 4^^^[cr3, i?^], so 
that 

Rj{x, ^, t) = e^'^'^'^Rjix, Oe-^'^'*''\ (20) 

Differentiate the equation (|T^) with respect to and then replacing by means 
of ( pO|) write it down in the form 

Applying the Liouville theorem to the function h and afterwards substituting 
the asymptotic representation ( p7| ) one can find that h takes the form Ai^^a^ — 
4:^^U — ^Vi — Vq (see formula (|T3D)- Consequently, solves the equation 

^j,t = 4^^'[^3, ^j] - {^eu + ^Vi + Vo)^j. (21) 

As it follows from equations (p^), ( |2T|) the function 

^(x,^,t) = ^,(x,^,t)e(4^^^-^^)-3 

is a solution to the over-determined system of equations (H^), (H^)? hence 
the potential U{x^t), found by solving the Riemann problem contains as a 
component a solution to the system of nonlinear equations (0), moreover, 
the initial conditions are automatically valid (|TT]). Check the validity of the 
boundary conditions ([TDD. To this end it is enough to show that at the point 
X = the system of equations ([T^) is invariant under rotation ^ cj^, i.e. 
not only the functions ^j{^,t), but also functions ^j{^,t) = ^j+2{<^^,t) are 
solutions to the system (PT|) . For values of the index j which are greater 
than six the function is defined periodically on j with period 6. By 
construction for any j the function $j(^, t) is analytic in the sector /j, as well 
as, the function mj(^, t) = ^'7^(^, t)^j{^, t) (requirements i)-iv), are supposed 
to still valid, therefore det^'j 7^ for ^ G Ij). As a result of this definition 
one gets a new Riemann problem similar (p^), 

^j+iiO = mVm, r,{0 = RM), (22) 
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with conjugation matrices rj{^) = mj{^)^^Rj{^)mj+i{^). When t = it can 
directly be shown that matrices mi, m3, are upper-triangular, and ma- 
trices 777-2, m4, 7776 are lower-triangular. By beans of this triangular structure 
the functions 

are also analytic in the corresponding sectors Ij. But it means that functions 
— ^) Q'lso solutions to the same Riemann problem 

(|22|), which is uniquely solvable, hence the functions ^j{^,t) and ^j{^,t) 
should coincide, consequently, m = m. Therefore the product {^jC^^^^^^^) x 
(^^.e4*^ ^3t)-i is equal to m(^,0), or in other words, it doesn't depend on 
i.e. functions y{^,t) = ^je'^^^^'^^^ and y{^,t) = ^je^^^^'^^^ solve one and the 
same equation of the form ([T^). By construction one has y{^,t) = y{u^,t), 
consequently, at the point a: = equation ( [T^ ) is not changed under rotation 
^ cj^, but then it is of the form (p7t|) . It is equivalent to the validity of the 
boundary condition (PD|). 

The conjugation matrices Ri and as triangular ones (see proposition 
2 and remark to it) remain analytical in the corresponding sectors I2 and 
for all t > 0. 

7 On solvability of the Riemann problem on a system 
of rays 

The Riemann problem on analytical conjugation on a system of rays which 
the inverse scattering is reduced to is not sufficiently well studied in spite 
of some close formulations has been met earlier in the connection with the 
Painleve equations But in our case the Riemann problem has inner 



symmetries which allow one to reduce it to the well known version of the 
problem posed on the real axis by transformations different for a: = and for 
X > 0. 

Let us begin with the case x = 0, t > 0. Note first that functions 
and ^5(^5 ^) are analytic not only in the corresponding sectors Ii and Iq but 
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also in larger domaims such as Ii U I2 and /s U /e, respectively, moreover 
they are related to each other by an involution ^5(^5 ^) = ^i{^^,t)- Define 
one more function by setting "^exi^^t) = ^i{^'^^,t)- These three functions 
constitute a solution to the following Riemann problem 

on a system of three rays ^i, ^3, £5 and in addition satisfy one and the same 
system of equations of the form Hence the conjugation matrices Mi, 

Mg, evaluate in time t as follows 

M(^, t) = e^'^'^'^Mi^, 0)6-^'^'"''. (24) 

It can be checked that the Riemann problems (|T^), ( PD|) and (|23|), ([24|) are 
equivalent: solution to one of them can easily be transformed into a solution 
to the other one by a simple explicit way. Let us given, for instance, a solution 
to the problem (|23|). Then functions ^3? ^4? ^5 which are the solution 
to the problem (|T5|) are expressed as 

where A^3 and A^4 are upper- and lower-triangular matrices. It is easy to check 
when t = that and A^4 are analytic in and I 4 and therefore, by means 
of the equations A^3^4(^,^) = e'^^^^'^^^N^^4{^,0)e~'^^^^'^^^ and triangular structure 
they preserve their analytical properties for t > also. 

Change variables ^ = i/^^^ in the problem (^). Evidently if ^ varies in the 
complex plane then u varies in three sheets Riemannien surface however in 
all these sheets the function defined as ^(i^, t) = ^i{^,t) takes one and the 
same values, i.e. it has no branch points, but fulfills a jump when crosses the 
positive half- line z/ > 0. Under the change of variables the Riemann problem 
on three rays is brought to the problem on real axis 

$(z/,t)*(z/,^) = C(z/,t), z/gR, (25) 

where ^{iy,t) = —A22{io^)'^Q^{^,t) and the conjugation matrix is chosen 
as C{iy,t) = -A22H) for z/ < 0, (^ G £2) C{iy,t) = -^2(c^O^6(0, ^, t) 
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for z/ > G £2)- Thus in order to solve the problem (|rB|) it is enough 
to solve the problem (|25|). Now it becomes clear that the problem (|T^) is 
solvable if functions 1/^2(0 ^ close to zero. This follows from the 

famous theorem about factorisation on a neighbourhood of the unity (see, 
for instance, [T^])- The positivity of the conjugation matrix also provides the 
unique solvability of the problem ([2^ ). Assume that the reduction constraint 
u'^{x^t) = —£u{x,t), £ = ±1 is imposed on the initial boundary value prob- 
lem (0), (pup, dTT]). Then columns i/j and (f) of the scattering matrix T satisfy 
the involution 02(0 = ^i(^*) = "^V^K^*)? where the star over a letter 

denotes the complex conjugation. In this case the conjugation matrix of the 
problem owing to the constraints i)-iv) from §5 is positively defined. 
Actually, the reduction condition gives Al2{^0 — ^22(^0 ^ ^ More- 
over, ^22(0) = ^22(00) = -1, ^22(0 hence, ^22(0 < along the ray £2 
such that C{iy,t) > for u < 0. Compute now the principal minors of the 
matrix C{iy,t): deti C{iy,t) = ^ii(0, det2 C{u,t) = det C{iy,t) = |^22(OP 
for 1/ > 0. By means of the constraint mentioned both are positive. Hence 
C{v^t) is positively defined and, besides, it is continuous in the Holder sense 
that provides the unique solvability of the problem ( |25|) . 

Consider now the case < x < 00^ t = const > 0, supposing that the 
Riemann problem (|T6|), ( pT|) , (p^ ), ( |2I]| ) is already solved for x = 0, t > 0. 
It means really that the scattering matrix T{^,t) = (?/^(^, t), 0(^, t)) of the 
problem on the half-line a: > is known for all ^ > and we come to 
the problem of recovering of the potential of the system ( [T^p for a given 
scattering matrix. The only solution of this problem (on the uniqueness see 
above proposition 1) can be found in two ways: by means of the Riemann 
problem ([T6]), (p7| ), ( p^) on a system of rays and by the canonical method 
based on the same problem on the real axis [T^: 

1 0i(^,i)exp(-2z^x) 
-^2(^,^)exp(2z^x) 1 
>,oo,^) = l. (26) 



(j)ix,^,t) = '0(x,<^,i) 



For a: > and for t > the Riemann problems mentioned are equivalent 
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unless the conditions i)-iv) are violated. Solutions to these two problems are 
related to each other by equations 

t) = t)mi{x, t), (j){x, t) = '03(x, t)m2,{x, t), 

^, t) = '05(x, t)mr,{x, t), t) = '(/^elx, t)mQ{x, t), 

where mj{x,^,t) are triangular matrix- valued piece-wise analytic functions. 
As an example we give an explicit representation of the matrix m^{x^ ^,t): 

m3(a:,^,t) = ^3V = 



,-12 1 ( Hi.t) 



It is well known (see, for instance, [1T2|) that under reduction u^(x,t) = 



—£u{x,t), £ = ±1 and the conditions i)-ii) the Riemann problem ([251) is 
uniquely solvable. 

It is worth to add that in the case e = —1 the conditions i)-iv) don't impose 
any constraint on the initial value of the initial boundary value problem (|T])- 
they are always valid. 
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